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Abstract I discuss constraints on star-formation and AGN activity in mas- 
sive galaxies at z~l-3 using observations from the Spitzer Space Telescope at 
3-24 /Ltm. In particular I focus on a sample of distant red galaxies (DRGs) with 
J — K s >2.3 in the southern Great Observatories Origins Deep Survey (GOODS- 
S) field. Based on their ACS (0.4-1 fxm), ISAAC (1-2.2 fj,m), and IRAC (3-8 fim) 
photometry, the DRGs have typical stellar masses A4>10 n A4q. Interestingly, the 
majority (>50%) of these objects have 24 /zm detections, with /„(24//m)>50 fiJy. 
If attributed to star formation, then this implies star-formation rates (SFRs) of 
| ~100-1000 Mq yr" 1 . Thus, massive galaxies at z~1.5-3 have specific SFRs equal 

to or exceeding the global average value at that epoch. In contrast, galaxies with 
A^IO 11 A4q at z^0.3-0.75 have specific SFRs less than the global average, and 
more than an order of magnitude lower than that at z~1.5-3. Thus, the bulk 
of assembly of massive galaxies is largely complete by z~1.5. At the same time, 
based on the X-ray luminosities and near-IR colors, as many as 25% of the mas- 
sive galaxies at z> 1.5 host AGN, implying that the growth of supermassive black 
holes coincides with massive-galaxy assembly. Lastly, the analysis of high-rcdshift 



galaxies depends on bolometric corrections between the observed Spitzer 24 /jm 
> ■ data and total IR luminosity. I review some of the sources of the (significant) 

uncertainties on these corrections, and discuss improvements for the future. 
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1 INTRODUCTION 

Most (~50%) of the stellar mass in galaxies today formed during the short time between z^3 
and 1 (e.g., Dickinson et al. 2003, Rudnick et al. 2003). Although much of this stellar mass 
density resides in massive galaxies, which appear at epochs prior to z~l-2 (e.g., Bell et al. 
2004, McCarthy 2004), it is still unclear when and where the stars in these galaxies formed. It 
may be that galaxies "downsize" , forming most of their stars in their current configuration at 
early cosmological times, with lower mass galaxies continuing to form stars to the present epoch 
(e.g., Bauer et al. 2005, Juneau et al. 2005). Alternatively, stars may form predominantly in 
low-mass galaxies at high redshifts, which then merge over time to form large, massive galaxies 
at more recent times (Bauer et al. 1998, Kauffmann & Chariot 1998, Cimatti et al. 2002). 

Understanding the formation and evolution of massive galaxies has been challenged by 
difficulties in constructing complete samples of star-forming and massive galaxies at z> 1. 
Massive galaxies at z<l exist on a fairly prominent red sequence (e.g., Blanton et al. 2003, Bell 
et al. 2004, Faber et al. 2005), and are largely devoid of ongoing star formation, evolve passively, 
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and contain up to half of the stellar-mass density. However, current hierarchical models predict 
colors for massive galaxies at z~0 that are too blue compared to observations (e.g., Somerville, 
Primack, & Faber 2001, Dave et al. 2005). Some recent models suppress star formation at 
late times in massive galaxies by truncating star formation at some mass threshold, or by using 
feedback from strong AGN (e.g., Granato et al. 2001, Dave et al. 2005, Croton et al. 2006). Most 
of the massive galaxies at z<l have red colors with formation epochs Zf> 2. The morphologies 
of the most optically luminous (and most massive) galaxies transforms from "normal" early- 
type galaxies at z~ 1 to irregular systems at z^ 2-3 (e.g., Papovich et al. 2005). Therefore, we 
need to study the properties of massive galaxies at these earlier epochs. 

In these proceedings, I discuss recent Spitzer observations at 3-24 ^im of massive galaxies 
at z~ 1.5-3 in the southern Great Observatories Origins Deep (GOODS-S) field. The Spitzer 
IRAC and MIPS observations provide constraints on the star-formation and AGN activity in 
massive galaxies at these epochs. I discuss implications the analyses of these data have for the 
stellar-mass assembly rates and formation epochs of massive galaxies. I also review some the 
uncertainties on these bolometric corrections between the observed Spitzer 24 ^im data and 
total IR luminosities in z^ 1.5-3 galaxies and prospects for future improvements. 

2 STELLAR MASSES AND STAR FORMATION IN HIGH -Z 
MASSIVE GALAXIES 

GOODS is a multiwavelength survey of two 10'xl5' fields, one in the northern Hubble Deep 
Field, and one in the southern Chandra Deep Field. The GOODS datasets include (along with 
other things) HST/ ACS and VLT/ISAAC imaging (Giavalisco et al. 2004), and recent Spitzer 
imaging. I make use of these data for the work described here, as well as data from Spitzer /MIPS 
24 ^m in this field from time allocated to the MIPS GTOs (e.g., Papovich et al. 2004). 

In these proceedings, I primarily focus on so-called distant red galaxies (DRGs) selected 
with J — K s >2.3 mag (see Franx et al. 2003). This color criterion identifies both galaxies at 
z~2-3.5 whose light is dominated by a passively evolving stellar population older than ~250 Myr 
(i.e., with a strong Balmer/4000 A break between the J and if s -bands), and also star-forming 
galaxies at these redshifts whose light is heavily reddened by dust (Forster-Schreiber et al. 2004, 
Labbe et al. 2005, Papovich et al. 2006). For the GOODS-S data, the J - K s >2.3 mag selection 
is approximately complete to stellar mass A4>10 n Mq for passively evolving galaxies, and we 
find 153 DRGs to K s <23.2 mag, spanning 0.8<z<3.7 with (z}~2.2 (see Papovich et al. 2006). 

More than 50% of the DRGs have 24 /im detections with /„(24/im)>50 [iJy. Daddi et 
al. (2005) find a similar 24 ^im-detection rate for massive galaxies at 1.5<2:<2.5 selected via 
their BzK colors. Interestingly, this implies that the majority of massive galaxies at z^2 emit 
strongly in the thermal IR — they are either actively forming stars, supermassive blackholes, 
or both at this epoch. The 24 /itm emission at z~1.5-3 probes the mid-IR (~5-10 /im), which 
broadly correlates with the total IR, Lir=L(8 — 1000/im) (Chary & Elbaz 2001). Figure 1 shows 
the inferred Lir for the DRGs using the Dale et al. (2002) models to convert the observed mid- 
IR to total IR luminosity. There is inherent uncertainty in this conversion, which I discuss in 
§ 3. The 24 ^m flux densities for the z-1.5-3 DRGs yield L m «10 n - 5 " 13 L Q , which if attributed 
to star-formation corresponds to SFRs of «100-1000 A4q yr _1 (e.g., Kennicutt 1998). 

Nearly all of the DRGs are detected in the deep Spitzer /IRAC data, implying they have 
substantial stellar masses. In Papovich et al. (2006), we modeled the DRG stellar populations 
by comparing their ACS, ISAAC, and IRAC [3.6yum] [4.5^m] photometry to a suite of stellar- 
population synthesis models (Bruzual & Chariot 2003), varying the age, star-formation history, 
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1 Total IR luminosities, 
1000), of galaxies inferred from their ob- 
served 24 fim emission. Red circles show 
the values for the DRGs; black stars show 
those objects with X-ray detections. Blue 
squares show objects from the HDF-N data 
(see Papovich et al. 2006). The solid line 
denotes the 24 /an 50% completeness limit. 
The inset error bar shows the estimated 
systematic error, «0.5 dex. 



Redshift 

Fig. 2 Stellar masses of galaxies inferred 
by fitting two-component models to the 
galaxies' rest-frame UV-to-near-IR data. 
The inset bars show the mean errors as 
a function of mass. The short-dashed line 
shows the characteristic present-day stel- 
lar mass (Cole et al. 2001); the long-dashed 
line shows the stellar mass limit for a pas- 
sively evolving stellar population formed at 
z^oo with K s —23.2 mag. 



and dust content. We use the model stellar-mass-to-light ratios to estimate the galaxies' stellar 
mass. We first allow for star-formation histories with the SFR parameterized as a decaying 
exponential with an e-folding time, r, ranging from short r's (burst of star-formation) to long 
t's (constant star-formation). We also use models with a two-component star-formation his- 
tory characterized by a passively evolving stellar population formed in a "burst" at Zf orm =oo, 
summed with the exponcntially-decaying-SFR model. Although the modeling loosely constrains 
the ages, dust content, and star-formation histories of the DRGs, it provides relatively robust 
estimates of the galaxies' stellar masses (see also Forster-Schreiber et al. 2004). Typical un- 
certainties for the stellar masses for the full DRG sample are 0.1-0.3 dex. Figure 2 shows the 
stellar masses inferred for the DRGs by fitting the two-component models. 

Figure 3 shows the specific SFRs / M) derived from the masses and SFRs for the DRGs, 
where the SFRs are calculated from the summed UV and IR emission. The figure also shows the 
specific SFRs for galaxies at lower redshift from COMBO-17 (Wolf et al. 2003), which overlaps 
with the GTO Spitzer 24 /im imaging. The SFRs for the COMBO-17 galaxies are calculated in 
an analogous manner as for the DRGs. The massive DRGs at 1.5<z<3 have high specific SFRs: 
DRGs with _M>10 n M e and 1.5<z<3 have */M~0.2-10 Gyr" 1 (excluding X-ray sources). 
In contrast, at z< 0.75 there is an apparent lack of galaxies with high specific SFRs and high 
stellar masses: galaxies with M>l0 n Mq have \E f /A / t~0.1— 1 Gyr -1 . 

We define the integrated specific SFR as the ratio of the sum of the SFRs, to the sum 
of their stellar masses, Mi, Y=J2i^i/J2i Mi, summed over all i galaxies. This is essentially 
just the ratio of the SFR density to the stellar mass density for a volume-limited sample of 
galaxies. Figure 4 shows the integrated specific SFRs for DRGs at z~l. 5-3.0 and COMBO-17 
at z ~ 0.4 and 0.7 with M > 10 11 M Q (see Papovich et al. 2006). The data point for the DRGs 
includes all objects with M>l0 n Mq, and assumes that 24 /im-undctectcd DRGs have no star 
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Fig. 3 The specific SFR as a function of galaxy stellar mass. The Left panel shows 
the DRG and HDF-N galaxies with 1.5 <z< 3, with symbols the same as in figure 1. 
The Right panel shows galaxies from COMBO-17 (as labeled). Open circles show 
COMBO-17 galaxies with 24 /im detections, small filled symbols show upper limits 
for galaxies undetected at 24 /iin. 

formation. The error-box lower bound shows what happens if we exclude objects with X-ray 
detections or IR colors indicative of AGN. The error-box upper bound shows what happens if 
we calculate SFRs for the 24 /im- undetected DRGs assuming they have /„(24/im)=60 pJy, the 
50% completeness limit. The integrated specific SFR in galaxies with Ai>10 11 M® declines by 
more than an order of magnitude from z^l.5-3 to z<0.7. Our results indicate that the relative 
star-formation in massive galaxies is reduced for 2<1 as galaxies with lower stellar masses have 
higher specific SFRs, supporting the so-called "downsizing" paradigm. 

Figure 4 also shows the specific SFR integrated over all galaxies (not just the most massive); 
this is the ratio of the cosmic SFR density to its integral, T=/3»/ J dt. There is a decrease in 
the global specific SFR with decreasing redshift. The evolution in the integrated specific SFR 
in massive galaxies is accelerated relative to the global value. Galaxies with .M^IO 11 Mq were 
forming stars at or slightly above the rate integrated over all galaxies at z~ 1.5-3. In contrast, by 
z < 1 galaxies with M. >10 n M.q have an integrated specific SFR much lower than the global 
value. Thus, by z<1.5 massive galaxies have formed most of their stellar mass, and lower-mass 
galaxies dominate the cosmic SFR density (see also Papovich et al. 2006). 

3 UNCERTAINTIES ON MID-IR DERIVED IR LUMINOSITIES 

3.1 AGN Contribution to the mid IR Emission 

Many of the massive galaxies at z~ 1.5-3 are detected in the deep X-ray data. In figures 1 
and 3, the X-ray-detected DRGs tend to have the highest inferred IR luminosities and specific 
SFRs. Many of these objects have Lir>10 13 Lq, comparable to PG quasars, which have warm 
thermal dust temperatures (Haas et al. 2003). The X-ray to optical flux ratios in these objects 
imply the presence of an AGN with Lx^lO 42 erg s^ 1 (for z >1.5, see Alexander et al. 2003, 
Papovich et al. 2006). In addition, many authors are finding AGN candidates based on rest- 
frame near-IR colors from Spitzer observations. The X-ray, UV, and optical emission in these 
AGN is heavily obscured by gas and dust, and they are often missed in deep X-ray surveys (e.g., 
Donley et al. 2005, Stern et al. 2005, Alonso-Herrero et al. 2006, Barmby et al. 2006). Roughly 
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Fig. 4 Evolution of the integrated specific SFR, the ratio of the total SFR to the 
total stellar mass (from Papovich et al. 2006). The curves show the expected evolution 
from the global SFR density (solid lines, Cole et al. 2001, thick line includes correction 
for dust extinction; dashed line, Hernquist & Springel 2003). Data points show results 
for galaxies with >10 n A4©. Filled circle corresponds to the DRGs; filled diamonds 
correspond to the COMBO-17 galaxies. The inset error bar shows an estimate on the 
systematics. 



25% of the non-X-ray detected DRGs in the GOODS-S field satisfy the IRAC color-selection 
for AGN from Stern et al. (2005). Of these, roughly one-half have ACS-through-IRAC colors 
consistent with dust-enshrouded AGN. Figure 5 shows examples of these galaxies with and 
with-out X-ray detections. In all cases, the galaxies show a clear flux excess in their IRAC 
[5.8/im] and [8.0/im] photometry, presumably arising from the AGN. Combined with the 15% 
of DRGs detected in the X-rays, up to 25% of the DRG population host AGN. 

If AGN contribute to the observed 24 /mi emission in galaxies at z^l.5-3, then they can 
affect the inferred IR luminosities. For example, although the Chary & Elbaz (2001) and Dale 
& Helou (2002) IR templates include galxaies with Ljr> 10 13 L®, using an IR template for 
Mrk 231 — with a known AGN and warmer dust temperature — would reduce the inferred IR 
luminosity for galaxies at 1.5-3 by factors of ~2-5. To limit the effects of any bias in the 
inferred IR luminosities caused by AGN, we restricted the error box on the integrated specific 
SFR for the DRGs in figure 4 to only those galaxies with no X-ray or IR indications for AGN. 

Even when AGN are present, it is unclear whether star-formation or the AGN activity 
domninates the bolometric IR luminosity. Alexander et al. (2005) demonstrated that a high 
fraction (~80%) of sub-mm galaxies are detected in 2 Msec Chandra X-ray data. However, 
the X-ray-detected sub-mm galaxies have IR to X-ray luminosity ratios up to an order of 
magnitude higher than what is expected for AGN alone. This suggests that both star-formation 
and AGN contribute to the bolometric emission. Similarly, Frayer et al. (2003) report a near- 
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Fig. 5 SEDs of DRGs with putative AGN. The data points in each panel show the 
rest-frame UV to near-IR SED from the ACS (0.4-l^m), ISAAC (1-2.2/xm), and 
IRAC (3-8/xm) photometry. The solid line shows the best-fit model SED to the ACS, 
ISAAC, and IRAC [3.6/xm] and [4.5/xm] photometry. The two galaxies in the Left 
panels have X-ray detections, but the two in the Right panels do not. In all cases the 
galaxies show excess emission in the IRAC [5.8/zm] and [8.0/im] bands, presumably 
from hot dust near the AGN. 

IR spectrum of a sub-mm galaxy SMM J04431+0210 at z~ 2.5. This galaxy has J - K s ~ 
3.2 mag, qualifying as a DRG. The spectrum shows that the galaxy nucleus has a low Ha to 
[Nil] flux ratio consistent with ionization from an AGN. However, Ha is spatially resolved in 
their spectrum and the [Nil] line strength drops off in the off-nucleus spectrum. This implies 
extended star-formation beyond the nucleus, which presumably contributes to the inferred IR 
luminosity. Therefore, it seems that both AGN and star-formation occur simultaneously in 
high-redshift IR-detected galaxies and both probably contribute to the IR emission. 

The high AGN occurrence in DRGs and sub-mm galaxies provides some evidence that 
massive galaxies at z~1.5-3 simultaneously form stars and grow supermassive black holes. 
Although speculative, the presence of AGN in these galaxies may provide the impetus for the 
present-day relation between black hole and bulge mass, and/or provide the feedback neces- 
sary to squelch star-formation (e.g., Kauffmann et al. 2004), moving the galaxies onto the red 
sequence. Emission line ratios from high-resolution spectroscopy at near-IR and/or mid-IR 
wavelengths may help constrain the ionization state of the galaxies' gas, identifying the fraction 
of galaxies with strong AGN activity as a function of mass at these redshifts. 

3.2 Uncertainties in the Shape of the IR Spectral Energy Distribution 

Although the mid-IR (5-15 /Ltm) emission broadly correlates with the total IR luminosity (e.g., 
Chary & Elbaz 2001), there exists considerable scatter. Partly this arises because galaxies of 
a given IR luminosity show a range in the strength of their mid-IR emission (from polycyclic 
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aromatic hydrocarbons, PAHs) and absorption features (e.g., Armus et al. 2004). This variation 
is reflected in a comparison of various model IR spectral energy distributions in the literature 
(e.g., see Le Floc'h et al. 2005). For example, if we used the IR templates of Chary & Elbaz 
(2001) instead of those of Dale & Helou (2002), then we would derive IR luminosities a factor 
of 2-3 higher for z^l.5-3 galaxies with Lir^10 125 ~ 13 Lq. In a recent study of galaxies at 
z<1.2 detected at ISO 15 /im and Spitzer 24 /<m in the northern GOODS field, Marcillac et al. 
(2006) found that, compared to Chary & Elbaz (2001), the IR models of Dale & Helou (2002) 
provide a tighter correlation between L\n derived from the mid-IR and radio flux densities with 
a scatter of 40%, suggesting the latter may better reflect reality. 

Some scatter is inherent in converting the 24 /im flux densities to total IR luminosity. 
Several studies of the mid-IR colors of galaxies to z<~l show that their 15 and 16 yum to 
24 (im colors have more scatter than predicted by simple models that map a single IR template 
to a given IR luminosity (e.g., Elbaz et al. 2005, Teplitz et al. 2005, Marcillac et al. 2006). 
Chapman et al. (2003) find that the temperature-luminosity distribution in local IR-luminous 
galaxies has a factor of 2-3 scatter in the IR-luminosity for galaxies of a given dust temperature. 
However, Daddi et al. (2005) found that the Chary & Elbaz (2001) IR model with L m =10 12 - 2 L Q 
fits the average spectral energy distribution of 24 ^m-detected BzK objects at (z) = 1.9 in 
the northern GOODS field. Therefore, while on average high-redshift galaxies have IR spectral 
energy distributions consistent with the models, individually there remains significant scatter. 



4 SUMMARY 

In this contribution, I discussed star-formation and AGN activity in massive galaxies 
(silO 11 -M©) a t z~l-3 using observations from the Spitzer Space Telescope at 3-24 /im. 
Interestingly, the majority (>50%) of these objects have /„(24/im)>50 /iJy, which if attributed 
to star formation implies SFRs of ~100-1000 M Q yr" 1 . Galaxies at z~l. 5-3 with _M>10 n M Q 
have specific SFRs equal to or exceeding the global average value. In contrast, galaxies at z~0.3- 
0.75 with these masses have specific SFRs less than the global average, and more than lOx lower 
than that at z~1.5-3. Therefore, by z<1.5 massive galaxies have formed most of their stellar 
mass, and lower-mass galaxies dominate the SFR density. 

As many as 25% of the massive galaxies at z>1.5 host AGN. The high AGN occurrence in 
massive galaxies at z^ 1.5-3 provides evidence that they are simultaneously forming stars and 
growing supermassive black holes. This may provide the impetus for the present-day black- 
hole-bulge-mass relation and/or provide the feedback necessary to squelch star-formation in 
such galaxies, moving them onto the red sequence. 

The largest source of uncertainty results from systematic errors on the bolometric correc- 
tions between the observed Spitzer 24 yum data and total IR luminosity. While on average 
high-redshift galaxies have IR spectral energy distributions consistent with the models, indi- 
vidually there remains significant scatter. Future work is needed to improve our understanding 
of the distribution between the mid-IR (rest-frame 5-15 /im) emission and total IR luminosity 
in galaxies at z~1.5-3. It may be possible to use "average" 24-to-70 and 160 yum colors of 
distant galaxies to further constrain the distribution and scatter of galaxies' IR spectral energy 
distributions. The upcoming Herschel Space Observatory (and eventually SAFIR) will mitigate 
this problem by measuring the far-IR emission of distant galaxies directly. 
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